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We introduce a systematic approach for the resummation of perturbative series which involve 
large logarithms not only due to large invariant mass ratios but large rapidities as well. Series of 
this form can appear in a variety of gauge theory observables. The formalism is utilized to calculate 
the jet broadening event shape in a systematic fashion to next to leading logarithmic order. An 
operator definition of the factorized cross section as well as a closed form of the next-to leading log 
cross section are presented. The result agrees with the data to within errors. 



Observables in weakly coupled gauge theories often ne- 
cessitate perturbative resummations to be under calcu- 
lational control. This need arises when one performs 
measurements that are sensitive to infrared scales. By 
probing distances long compared to the hard scattering 
scale one introduces large logarithms (logs) that lead to 
the breakdown of fixed order perturbative series. Re- 
summing the large logs has become standard in QCD [1 
and can be accomplished by factorizing the cross section 
into momentum regions. Factorization makes clear the 
distinction between logs of various ratios that may be 
involved in the observable and resummation follows via 
standard renormalization group techniques. 

An elegant formalism for factorization is SCET (Soft- 
Collinear Effective Theory) [2] , which is an effective field 
theory designed to reproduce the infra-red physics in high 
energy processes. The formalism not only streamlines 
factorization proofs [5] , but also allows one to systemati- 
cally include power corrections. The results of this paper 
will all be couched in terms of this framework. 

A generic factorized cross section takes on the form: 

a = H®[Il i J i ]®S. (1) 

The hard function H is responsible for reproducing the 
short distance physics with wavelengths of order 1/Q, 
where Q is the scale involved in the hard scattering. J.; 
and S are the so-called jet and soft functions containing 
modes which are highly energetic (collinear) and soft, re- 
spectively. Soft modes have small rapidities (fc+/fc_ ~ 1) 
while the rapidities of collinear modes are parametrically 
larger (k±/k T 3> 1), where k± are the light cone mo- 
menta. The tensor product implies the existence of one 
or more convolution in momentum space. In canonical 
situations, the resummation of large logs is accomplished 
by evolving, via the renormalization group (RG), each 
factorized component to its natural scale. The natural 
scales are set by the arguments of the logs. H, J and S 
may contain, for example, logs of Q/ n,mj/ )u,ms/7x, re- 
spectively, where mj and ms are quantities which probe 
the invariant masses of the modes composing J and S. 

While there is a large disparity in rapidity between the 
modes which compose S and J, the typical invariant mass 
of the modes need not be distinct. When soft fields have 
invariant mass parametrically smaller than the collinear 



modes (in this case the soft mode are called "ultra-soft" ) 
significant simplifications arise. Whether or not there is 
a distinction in invariant masses, one must always ensure 
that there is no double counting of modes. That is, loop 
integrals within a prescribed function (J or S) should 
only account for the relevant mode. In principle this 
could be accomplished using a cut-off, but this would lead 
to multiple technical difficulties not the least of which is 
the need for gauge non-invariant counter-terms. Within 
the effective field theory formalism using dimensional reg- 
ularization, this double counting is avoided by the so- 
called zero-bin procedure [I]. In this methodology one 
subtracts from each loop integral its value when the inte- 
grand is asymptotically expanded around the extraneous 
region. This procedure not only formally avoids the dou- 
ble counting, but also ensures that all integrals in the 
theory are well defined. Moreover, the zero bin subtrac- 
tion, or some equivalent subtraction method, is necessary 
to preserve factorization 5J. This potential breakdown 
of factorization occurs as a consequence of the need to 
regulate "rapidity divergences" (light-cone singularities) . 
These divergences arise schematically from integrals of 
the form 

which are not regulated in dimensional regularization. 
There are multiple ways of regulating this divergence. A 
simple way is to introduce a new dimensionful parameter 
A via the replacement k + — > k + + A in the denominator. 
Note that regulating these divergences will break boost 
invariance along light-cone direction. For any physical 
observable the final result must be boost invariant and 
independent of A, this is automatic once zero-bin sub- 
traction is performed and all sectors added. For the case 
of to j 3> rn s , A dependence cancels in each sector after 
zero bin subtraction, and the boost symmetry is restored 
in each sector. This will, however, not be the case when 
to j ~ ms, since the soft and collinear modes mix under 
boosts. 

When toj ~ to 5 the jet axis, defined as the direction 
of the net momentum of the jet, recoils against the soft 
emission. In the light-cone co-ordinates, collinear modes 
scale like (n • p,n ■ p,p±) ~ Q(l, A 2 , A) while soft modes 
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FIG. 1. Rapidity renormalization group flow along the on- 
shell hyperbola versus the standard flow. 



scale as Q(X, A, A), where A< 1 and is the light-cone 
direction chosen to perform the factorization. The jet 
axis is no longer aligned with n M , and one should not ex- 
pect the jet function to be invariant under boosts along 
. However, the sum of all sectors will still be invariant. 
If we were to reanalyze the rapidity divergences discussed 
previously we would find that the A regulator will cancel 
after summing over contributions from all sectors (with 
proper zero-bin subtractions) [5 ]. Boost non-invariance 
in the jet function appears in the form log[Q/A], a "ra- 
pidity log" . Resumming these logs using a A regulator, 
is technically cumbersome much like resummations with 
a cut-off regulator. Here we introduce a regulator more 
in the spirit of dimensional regularization that does not 
introduce new dimensionful scales in the integrals, and 
maintains manifest power counting in the effective the- 
ory. 

Given the existence of the rapidity logs in addition to 
the canonical logs, S and J may not have one definitive 
scale associated with them. To resum both sets of logs we 
will introduce another arbitrary scale i/, along the lines of 
\x in dimensional regularization. We expect that in order 
to properly resum all the large logs we will need to run 
the jet in v down to the smaller rapidity scale of the soft 
function. In a Wilsonian sense we have two distinct cut- 
offs with which to thin degrees of freedom. There will be 
one flow in invariant mass and one in rapidity as shown 
in Fig. [l] This is an inherently Minkowskian procedure. 

To illustrate this rapidity renormalization group we 
will consider the specific example of the event shape 
called jet broadening. Event shapes have played an im- 
portant role in precision measurements of the strong 
coupling a s [6 . A generalized event shape for event 
e~e + — > X at center of mass energy y/s, can be defined 
[7] in terms of a parameter a via 



e(a) 
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u J±i e -\m\{i- 



(3) 



where Pi± is the transverse momentum with respect 
to the thrust axis of the event, and r\i is the rapidity 
of the i'th particle. The thrust axis t is defined via 
max t - X^iex \Vi'AI \f^- Two particularly interesting event 
shapes are the limits a = 0, 1 corresponding to "thrust" 
and "jet broadening" respectively. The limit e -C 1 iso- 
lates events composed of back to back jets. In the case of 
thrust these jets are composed of collinear radiation, and 
the recoil due to soft (ultra-soft in this case) radiation 
does not affect the jet axis. While for jet broadening all 
radiation with parametrically similar transverse momen- 
tum can contribute, so that soft radiation of the form 
Q(X, A, A) recoils the jet off the thrust axis. In both of 
these cases fixed order perturbation theory will fail when 
e is small. However, as long as eQ Aqcd, we expect 
non-perturbative effects to be suppressed, though large 
logs of e need to be resummed. 

The pioneering work on jet broadening resummations 
[8] utilized the coherent branching formalism [10] . It was 
later stated [5] that the results in [5j neglected terms due 
to recoil of soft gluons. In this letter we will provide 
a factorization theorem for jet broadening, whose proof 
will follow in a subsequent publication [H]. The factor- 
ization proofs for angularity observables ([3| in pj] are 
known to fail as a approaches 1, since there are growing 
power corrections in this limit where one approaches jet 
broadening. The reason for the apparent breakdown of 
factorization is the fact that in this limit the soft radi- 
ation has the same invariant mass as collinear radiation 
and one must change the power counting accordingly to 
factorize in a consistent fashion [TJ]. 

Hereforth, we set a = 1 and e(l) = e. In [T3] we 
prove that the cross section for jet broadening takes the 
following factorized form 



da 
de 



= a Q H(s) / de n de n de s 8(e - e„ - e n - e s 



(4) 



J dp lt dp2tJn{Q+,e n ,Plt)Jn(Q-,e rl ,,P2t)S(e s ,p 1 t,P2t) 



where in covariant gauges 
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Xn(x + y-S(e-e n )5(P x +p lx )xn(0)\0), 



s = p\^ P ^n d f 



(5) 



(0| SiS n 5{e - e s )S d (P nl _ - pii_)6 d (P nl _ - p 2 ±)Sts n \ 0), 

and (To is the Born cross section. H is the hard matching 
coefficient which incorporates all the short distance con- 
tributions, and is fixed by matching the QCD currents 
onto the SCET currents. Here, d — 2 — 2e, X(n,n) are 
gauge invariant SCET fields which include collinear Wil- 
son lines W n> n and the S n ^i are soft wilson lines. P n ±_ and 
PfiX are hemisphere-transverse momentum operators, fig 
refers to area of a (i-dimensional unit sphere, fl±2 the rela- 
tive angle between p%± and p2±_ , and pu = \pn__ \ ■ Finally, 
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Q ± are the large light cone momenta of the jets with 
constraint Q + Q~ = s, the center of mass energy. 

As long as y/se ^> Kqc d all of these matrix elements 
are calculable in perturbation theory. The bare matrix el- 
ements possess both rapidity and UV divergences. Thus 
standard dimensional regularization is insufficient to reg- 
ulate all the integrals. Beyond tree level one is immedi- 
ately met with the aforementioned rapidity divergences. 
To regulate these integrals we introduce a regulator into 
the momentum space SCET Wilson lines as follows 



W n = 



J2 ex p 



\n-Pr v 
w l - l —n-A n , q {0) 



perm 



E 

perm 



exp 




n-A s , q (0) , (6) 



where v is an arbitrary scale independent of the usual 
scale [i introduced in dimensional regularization. to is a 
bookkeeping parameter which will be set to one at the 
end of the calculation. P M is the momentum operator 
and we have essentially regulated the energy of the emit- 
ted gluons in each Wilson line, since 2Po — > n. • P in the 
collinear limit. Notice the factor of r//2 in the soft func- 
tion. This choice is not a matter of convention. Physi- 
cally, the factor arises as a consequence of the fact that 
soft must be cut-off at both positive as well as negative 
rapidity [14] . The Wilson line regularization breaks man- 
ifest boost symmetry, which is restored once all of the 
sectors are combined. The rapidity divergences for the 
jet and soft functions will introduce a new set of anoma- 
lous dimensions (7^, 7^) which are defined via variation of 
v. Given that the hard function has no such anomalous 
dimensions, we must have the relation 



27^ + 7s = 0, 



(7) 



or equivalently it must be true that the total r\ depen- 
dence must vanish. Indeed it is not hard to show |14j . 
that the sum of the 77 divergences cancels as a conse- 
quence of eikonalization. This cancellation also implies 
that the individual factors J and S are multiplicatively 
rcnormalizable . 

The tree level jet function is given by 5{e n — pu/y/s) 
and soft function by S(e s )S(pit)S(p2t)- To determine the 
relevant scales in the logs, we must convolve the renor- 
malized one-loop jet function with the tree level soft func- 
tion as dictated in Q. This is most easily seen for the 
integrated cross-section S = J Q e ° de (da / de) . At order 
a s , the result of the convolution leads to the following 
singular contributions from a jet 
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It thus becomes clear that the jet function depends on 
two different kinematical scales (\/seo and Q ). In ad- 
dition, we see that dependence on rapidity manifests in 
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FIG. 2. Running Strategy. 



the form of the rapidity log, \m(Q ± /v). While the soft 
function singular contributions are 
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and has rapidity logs set at the low scale y/seQ. 

The utility of v is clear as we may choose v ~ Q ± ~ \fs 
and /i ~ yfse to minimize the logs in the jet function. 
Then to minimize the logs in the soft function we run v 
from the scale y/se up to y/s. Furthermore we will need 
to run the hard matching coefficient down to the scale 
\J~se. in the dimensional regularization parameter \i. This 
scenario is shown schematically in Fig. [2j 

Here we will perform the running at next to leading 
log (NLL) , which sums all terms of order one, where 
we take the scaling a s ln(e) ~ 1. The NNLL analysis 
will be performed in [TJ. The hard function RGE is 
well known to NNLL (see [Hj). The running is most 
simply performed in Laplace transform space, with (6, b') 
conjugate to (puiPit) respectively. We find 

6') = -2r c [(log(&/^ B ) + log(&V 7B )] , (10) 

where r c is the usual cusp anomalous dimension for Wil- 
son lines. Similar equations can be written for the two 
jet functions in terms of their corresponding anomalous 
dimensions. 

In our strategy of resumming logs for the soft function 
we only need two-loop cusp. Solution of the I/-RGE for 
the soft function is 



S(fi, v) = V s (n, v/v ) ® u ) 



(11) 



where ® represents convolution in kinematical arguments 
which are dropped for brevity. Here, 



Vs(pit,P2uU s ([J,,ls/vo)) 
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(12) 

S(P2t) ■ 



with u> s (fj,,v/uo) = 2T c [a s (p,)] log ^ . To minimize logs 
in the hard function we need to evolve the hard function 
using 
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H(s; fx) = H(s; no) U H {s; juO, fj.) 



(13) 



4 



130 GeV 




® Data from L3 

::::: LL 

NLO sing. 



0.1 0.2 0.3 Bj 

FIG. 3. Total Jet Broadening at 130 GeV. 



where up to NLL Uh can be found in |12j . and H(s, Ho) = 
1 to the order we are working. 

The results we have presented so far are for the angu- 
larity at a = 1, which is related to the total jet broad- 
ening Bt via e = 2Bt- We will present cross-sections 
for total jet broadening here and compare with the data. 
For NLO singular cross-section we get 



da 
~d~Eh 



a s (ij)C F 

00 5 (- 

7r Bt 



3-41og J B T ) 



(14) 



where cxo is the Born cross-section. This result is in agree- 
ment with ref. [5]. For the re-summed cross-section up 
to NLL order we have 



da _ a U H {Q 2 , HQ, fJ.) ( QB T 



dB T B T T{2u s )e 2 "f^° 



(15) 



In Fig. [3] we have plotted the theory cross section and 
the data [13]. We see that given the large error bars 
the agreement with the data is reasonable. However, the 
NNLL calculation will reduce the theory errors consid- 
erably. We have not included the theory errors due to 
power corrections. In the small Bt region these are non- 
perturbative and scale as ^-qcd/(BtQ) and can be ex- 
pected to be of order 20-30%. In the tail region there 
are corrections of order Bt relative to the singular con- 
tributions. The disagreement at intermediate values of 
Bt, where fixed order calculations suffice, is expected, 
since logs will not dominate in this region and NLL re- 
sults leave off order one contributions. This region will 
be correctly reproduced in the NNLL calculation. There- 
fore, by systematically improving this result by includ- 
ing higher order corrections in a s , power corrections and 
non-perturbative correction, this result can be used for 



precision a s determination. Such an analysis using thrust 
was done in [J"2] . 

Finally, we wish to point out that the rapidity renor- 
malization group can be utilized in multiple other set- 
tings where rapidity divergences arise. Generically, this 
will occur whenever kinematically soft radiation has in- 
variant mass of the same order as the collinear radiation, 
as in cases where one measures the pt of the final state. 
Such observables will be discussed in more detail in |14j . 
Furthermore, it would be interesting to utilize our ra- 
pidity renormalization group in the context of exclusive 
processes where it has been shown rapidity factorization 
sheds light on end point singularities in integrals over 
light-cone distribution functions [4]. 
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